Tunneling of quantum rotobreathers 
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We analyze the quantum properties of a system consisting of two nonlinearly coupled pendula. 
This non-integrable system exhibits two different symmetries: a permutational symmetry (permu- 
tation of the pendula) and another one related to the reversal of the total momentum of the system. 
Each of these symmetries is responsible for the existence of two kinds of quasi-degenerated states. 
At sufficiently high energy, pairs of symmetry-related states glue together to form quadruplets. We 
show that, starting from the anti-continuous limit, particular quadruplets allow us to construct 
quantum states whose properties are very similar to those of classical rotobreathers. By diagonaliz- 
ing numerically the quantum Hamiltonian, we investigate their properties and show that such states 
are able to store the main part of the total energy on one of the pendula. Contrary to the classical 
situation, the coupling between pendula necessarily introduces a periodic exchange of energy be- 
tween them with a frequency which is proportional to the energy splitting between quasi-degenerated 
states related to the permutation symmetry. This splitting may remain very small as the coupling 
strength increases and is a decreasing function of the pair energy. The energy may be therefore 
stored in one pendulum during a time period very long as compared to the inverse of the internal 
rotobreather frequency. 
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I. INTRODUCTION 



As revealed by the increasing number of recent pa- 
pers [|l|-[^6|, the quantum counterpart or the quantiza- 
tion of classical discrete breathers (CDBs) has become 
these very last years a real challenging and exciting field. 
At present, indeed, the theory of CDBs, i.e. "time- 
periodic spatially localized motions in networks of oscilla- 
tors (MacKay |2j|)" , has reached a high degree of perfec- 
tion and may certainly be considered as a real achieve- 
ment. As it is well known, the necessary condition re- 
quired for a Hamiltonian lattice to support such CDBs 
is that the breather frequency as well as its harmonics 
do not lie inside the phonon band. Such a non resonance 
condition may be achieved by the interplay of discrete- 
ness and nonlinearity, the first providing a natural up- 
per bound of the linearized Hamiltonian spectrum (say 
phonons) whereas the latter allowing for a tuning of fre- 
quencies out of this spectrum Q . 

Mathematical existence proofs of CDBs have been ob- 
tained for a wide variety of model Hamiltonians pq ] 
and have simultaneously given rise to an accurate nu- 
merical method allowing for their practical construction 
p9| . Both are based on the so-called anti- continuous 
limit (i.e. the limit where the coupling between oscil- 
lators vanishes), provided a nonlinear onsite potential is 
present. In this limit time-periodic spatially localized (on 
one or more sites) solutions trivially exist. More recently, 
an existence proof of CDBs based on a discrete version 
of the center manifold reduction has also been given for 
the FPU chain for which the method based on the anti- 
continuous limit fails [Q. This completes and enlarges 
the results previously obtained in |31| by a homoclinic 
orbit approach. At last, also recently, a variational ap- 
proach has been carried out to prove the existence of hard 



discrete breathers in some classes of Hamiltonians [ p2[ . 
Obviously, the wide range of applications of these rigor- 
ous mathematical results supports the idea that CDBs 
are generic solutions of nonlinear Hamiltonian lattices as 
claimed for the first time by Sievers and Takeno in 1988 
when they discovered this new kind of intrinsic localized 
modes [ p3| . 

Together with their existence proofs, the properties of 
CDBs have been extensively studied. It has been shown 
for example that they are structurally stable provided the 
non resonance condition holds ]27j and linearly stable in 
any dimension provided the coupling is weak enough p[ . 
Their spatial decay is generally exponential (for any finite 
range interaction potential) but may be also algebraic in 
case the interaction potential itself is algebraically de- 
caying p3| . More recently, it has even been confirmed 
that a pure nonlinearity of the interaction potential may 
give rise to a super-exponential decay of the breather 
tail |27], |3f|. At last, and directly related to their spatial 
decay properties, possible energy thresholds for their ap- 
pearance have been derived according to whether the lat- 
tice dimension exceeds a system-dependent critical value 
or not |36| . Extensive studies of CDBs properties as well 
as the perspectives in this field may be found in different 
reviews (see for instance Q, [l5|, |2j| |2j| |37j ) . 

The general properties of CDBs briefly listed above 
shows how well their theory is understood at present. 
Much less is known on the other hand about their quan- 
tum counterparts (QDBs), that is, about the quantum 
states which would behave similarly to CDBs in a sense 
to be clarified below. Nevertheless, this question is of 
great importance given one natural area of application, 
namely, the condensed matter physics where quantum ef- 
fects have generally to be considered (see |23] and refs. 
therein for a review of the materials where the presence of 
DBs is suspected, expected or has been detected). This is 
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the precise reason why the quantum theory of breathers 
is currently in full development. 

Let us turn back to an essential property of CDBs. 
These time-periodic solutions occur in a lattice but are 
spatially localized which means that they are not invari- 
ant under the discrete translational symmetry of the lat- 
tice whereas the Hamiltonian does. If we now consider 
the corresponding quantum problem, the invariance of 
the quantum Hamiltonian H with respect to the discrete 
translational symmetry (the operator of which we de- 
note by T) yields [H,f] = 0. The eigenstates of f are 
thus eigenstates of H and are delocalized along the lattice 
(Bloch waves). Of course, this discrete translational sym- 
metry of the lattice is broken if we consider for instance 
a finite system with fixed or open boundary conditions. 
Nevertheless, if the lattice consists of a sufficiently large 
number of sites, discrete translational symmetry is prac- 
tically restored similar to the infinite lattice. However, 
at large enough energy the quantum and the classical 
descriptions of the system should give similar results! 

It is possible to reconcile the quantum and the clas- 
sical points of view provided we again start from the 
anti-continuous limit. In this limit indeed, the Hamil- 
tonian consists only of a sum of identical Hamiltonians 
(one for each site) H = J2 S Hs and a general eigenstate 
\t/j) of the system may be represented as the tensorial 
product of local eigenstates |0 S ), \ip) = ® s \4> s ). Then, 
a localized excitation is obviously constructed by cre- 
ating an excitation of level k at site n, \4>^} whereas 
the other sites reside in their ground state \4>f^ n )- This 

leads to = ® s <n\^s ] ) <8> \<t> { n) ®s>n \<f>^). This 

case corrresponds to the classical one used to initialize 
the construction algorithm of a CDB at non zero cou- 
pling. Nevertheless, any eigenstate of the form ) — 

(3s<]\<f>^} <8 \<t>^) ®s>j |<^) has the same energy. The 
corresponding eigen-subspace has a dimension N equal 
to the number of sites (provided no additional accidental 
degeneracy occurs). Therefore in the uncoupled (anti- 
continuous) limit, any unitary transform of the proceed- 
ing basis leaves the subspace invariant although it yields 
a new basis. As the Hamiltonian system is invariant un- 
der the discrete translation T whatever the intersite cou- 
pling e is, at non zero coupling, the eigenstates of H must 
belong to one of the N symmetry sectors defined by T. 
When the coupling becomes zero, we may thus choose 
the eigenbasis of T as a basis of the ./V-dimensional sub- 
space. We denote its eigenvectors by |Ag fc ^) where q la- 
bels the symmetry sector. Of course each of these new 
eigenstates is completely delocalized (in the sense that 

v<z, E.KVf lAf)! 4 ]- 1 ^)- 

Let us assume now that we switch on the coupling 
between the sites. We expect the iV-fold degenerated 
eigenenergy to be splitted under the effect of the per- 
turbation and to form a band of N nearly degenerated 
eigenenergies. The corresponding eigenstates, which al- 
ready had the correct symmetry because of the unitary 



transform T performed on the local basis at e = 0, now 
give a new (perturbed) basis \Xq k \e)). Provided e is 
weak enough, these new eigenvectors are close to those 
defined at the uncoupled limit and thus an inverse uni- 
tary transform T _1 is expected to yield a basis \ip^\s)) 
close to the local basis of the uncoupled limit, that is, 
|^ fc) (e)) = \ip ( n k) )+0(e). These states are thus well local- 
ized provided the coupling is weak enough. However, be- 
ing a linear combination of non degenerated eigenstates, 
they are no more eigenstates of H and evolve in time. 
This tunneling effect, corresponding to the transfer of an 
initial excitation from site to site, takes a time typically 
given by the width of the band at a given value of the cou- 
pling s. In order to gap the bridge between the classical 
manifestation of the breather solution and its quantum 
realization, we thus expect the quantum system to ex- 
hibit bands of N nearly degenerated eigenstates whose 
width tends to zero as the average energy of the band 
goes to infinity. 

The purpose of a quantum theory of breathers is thus 
to know how the above defined bands behave as the cou- 
pling increases. Several successful attempts to answer 
this question have already been done in one-dimensional 
systems which present the particularity to be integrable 
(@, H-0 and the review §|). The (soliton/breather) 
bandwidth is shown to behave typically like AE n ~ 
(e/j) n /(n — 1)! where n is strictly speaking the num- 
ber of bosons in the system. This number corresponds 
to the number of the excited level of a single site in the 
absence of interaction e — and typically depends on en- 
ergy as a power law. Finally 7 is a measure of the onsite 
(quartic) nonlinearity. It is then obvious that for bands 
of high energy (n — > 00), the bandwidth goes rapidly (in 
fact more than exponentially fast) to zero with increasing 
n. In this case the tunneling time goes to infinity in the 
classical limit and the quantum breather (quasi-classical) 
state remains localized on its initial site. 

Several questions remain however. At high energy, if 
the onsite potential behaves like x 2q , the density of states 
(evaluated at e = by the mean of the Weyl's formula) 
scales like g(E) ~ a s E -> 00, where N 

is the number of sites. An immediate conclusion is that 
if the number of sites is greater or equal to 2, the den- 
sity of states (DOS) increases with energy according to 
a power law with exponents linearly depending on N . 
Then, we may ask whether the corresponding raise of 
interaction possibilities (hybridization) will destroy the 
quantum breather bands or not. A partial answer to 
this question should be given by the ratio of the breather 
bandwidth and the mean level spacing l/g(E) at a given 
(high) energy E: lim^^oo g(E) AE ? Assuming that the 
DOS scales as indicated by the Weyl's formula, the ra- 
tio of the bandwidth to the mean level spacing will tend 
to zero provided the bandwidth decreases exponentially 
with increasing energy (or faster). The above mentioned 
integrable examples fall into this category. This result 
seems to indicate and to explain the possibility, starting 
from the quantum problem, to recover the limit of the 
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classical breather at high energy. Another natural ques- 
tion is related to the impact of the nonintegrability of 
the system on the quantum breather bands and in par- 
ticular the influence at the quantum level of the chaotic 
trajectories induced by the nonintegrability lying nearby 
the classical breathers. 

Up to now, except for integrable systems, where 
analytical results are obtainable, studies of quantum 
breathers have been done numerically. Such studies be- 
come rapidly extremely difficult due to the huge matri- 
ces to be diagonalized and so far have been restricted to 
small ID systems for which the number of sites has not 
exceeded N = 12 ^0|. Even for these moderate lattice 
sizes, the average dimension of the (truncated) Hamilto- 
nian matrices is generally of order 10 6 and req uires spe- 
cific numerical diagonalization methods Jl3| , [lq |. The nu- 
merical results reported in the papers mentioned above 
are restricted to the low energy sector as soon as the 
number of sites exceeds a few units. This is easily un- 
derstood given the number of configurations obtainable 
by truncating the basis to p bosons per sites on N sites 
((p + 1) ). Remains the possibility to study very small 
systems (N — 2 or 3), as it has already been done for 
the dimer (integrable) [Q, [l2| and the trimer (non inte- 
grable) |55| . The dimer has also been used to de- 
scribe the tunneling of a QB along a chain by a suitable 
linearization of the lattice around it 18 . This lineariza- 
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tion method was also employed to study the properties 
of classical rotobreathers in a chain of pendula [^2| . But 
the specific properties of quantum rotobreathers have not 
been discussed yet. We will consider a dimer of two cou- 
pled quantum pendula which, due to the presence of both 
nonlinear onsite and interaction potentials, is not inte- 
grable. 

As far as we know, no experimental study of quan- 
tum rotobreathers has been done yet. This is possibly 
due to their quite large activation energy. The reader 
should notice however that a large number of studies have 
been devoted to the rotational motions of molecules (see 
e.g. @, 



especially concerning the methyl 
groups whose dynamicalproperties are usually obtained 
via neutron scattering |f § H, |§. But the tran- 
sitions thus obtained concern the so-called quantum ro- 
tational tunneling effect (43). This tunneling occurs be- 
tween the equilibrium positions defined by the onsite po- 
tential which is ?i-fold according to the individual sym- 
metry of the observed molecule and the symmetry of its 
environment. For the methyl groups n is generally equal 
to 3 but for ammonia in Hofmann clathrates, the four- 
fold symmetry of the host crystal induces an approxi- 
mative 12-fold |44| or more complex [Ell ^] symmetry. 
This tunneling is thus responsible for a rotation of the 
molecule forbidden in the classical case. The quantum ro- 
tobreather on an other hand is a state whose energy is sit- 
uated above the energy of the separatrix. A tunneling ef- 
fect appears, because of the coupling between molecules, 
which consists in the transfer of the excitation from site 
(molecule) to site. A study of the properties of the 4- 



methyl-pyridine by Fillaux and co-workers has revealed 
the presence of a quantum sine-Gordon breather in this 
compound M, [|, [l7j . The properties of this state, shown 
to be the ground state of the system, are theoretically an- 
alyzed via a semiclassical quantization procedure of the 
classical solution of the sine-Gordon equation and then 
successfully compared to experimental results. However, 
the coupling between adjacent methyl groups is so strong 
along certain axes of the crystal (chains) that the relative 
phases on neighbouring groups are small. That allows to 
use the continuum sine-Gordon theory. In the case of a 
rotobreather, such an approximation becomes of course 
impossible due to the unavoidable large phase difference 
created at the interface between the rotating and oscillat- 
ing groups. The study of the properties of rotobreathers 
(either classical or quantum) thus requires to preserve 
the rotational invariance of the interaction potential. 

Our system under study consists of two pendula cou- 
pled by a cosine interaction potential (which preserves 
the rotational invariance). To simplify the study, the 
onsite and interaction potentials are the same functions 
and are 1-fold. This avoids to deal with the nearly degen- 
eracies coming from the n different arches of a cos(nx) 
potential at low energies (typically below the separatrix) . 
Notice that this does not change essentially the high en- 
ergy part of the single site spectrum (which tends to the 
free rotor one whatever the value of n) but removes the 
splittings of the doublets which are not coming from the 
center (k = 0) or the edge (k = n/2) of the corresponding 
Brillouin zone k 6 {— n/2 + l,..,n/2}. The results we will 
obtain in section 2 concerning the upper part (above the 
separatrix energy) of the spectrum are thus easy to adapt 
to the case of an n-fold onsite potential. Notice at last 
that despite its very short size, our "pendula dimer" still 
possesses the discrete translational invariance T which 
corresponds to the permutation symmetry. 

The plan of the paper is as follows: in section 2 we 
briefly review the properties of a single quantum pendu- 
lum and we derive important formulas concerning the 
splitting occurring in the doublets in the high energy 
sector. In section 3, we start by presenting some re- 
sults of the classical problem corresponding to two cou- 
pled pendula. Then we show how to compute the quan- 
tum spectrum of this system. We derive the exact spec- 
trum at the uncoupled limit and we construct a quan- 
tum rotobreather according to the method presented in 
the introduction. By increasing the coupling between the 
pendula, we follow the quantum rotobreather (QR) and 
we compute the splittings occurring in the corresponding 
quadruplet. 



II. SINGLE PENDULUM PROBLEM (SPP) 

In this section we briefly review some of the proper- 
ties of the classical and quantum one pendulum problem 
(for more details see for instance [p6j| ) . In the quantum 
case, we particularly focus on quasi-degenerated states 
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and their energy splitting. We show by two different 
methods that it is possible to compute this splitting ex- 
actly in leading order. The Hamiltonian of the pendulum 
system is given by 



H(x,p) — \-p 2 + a(l — cos(x)) 



(1) 



where x and p represent respectively the angle variable 
and the associated momentum, a > tunes the barrier 
height of the onsite potential. 



A. Classical case 

From ([!]) we deduce the classical equation of motion 

x + asin(x) = (2) 

which possesses two different solutions. Denoting by E 
the energy of the pendulum and by E s — 2a the energy 
of the separatrix separating the oscillatory motion from 
the rotational one we get: 



• Oscillation: E < E s 

sin(a;/2) = fcsn {yfctt\ k) 



(3) 



where the modulus of the Jacobian elliptic function 
sn is defined by the relation k = (E/2a) 1 ^ 2 and the 
period of the oscillation is T osc = AK(k)/ \fa 1 K(k) 
being the complete elliptic integral of the first kind 



• Rotation: E > E„ 



sin(a;/2) = sn ( —~-t;k 



(4) 



where k = (2a/ E) 1 / 2 and where the rotation period 
is T rot = 2kK(k)/^a~. 

B. Quantum case 

Surprisingly, while the classical problem of the pendu- 
lum is a very basic one, it seems that a thorough study 
of the related quantum problem has been done quite re- 
cently (R. Aldrovandi and P. Leal Ferreira fi(|). Our 
purpose here is not to repeat the analysis done in the 
paper cited above but rather to insist on the momentum 
reversal symmetry of the Hamiltonian (|l]) which leads to 
the appearance of pairs of quasi-degenerated states above 
the separatrix. 

1. Analytical solution 

The stationary Schrodinger equation corresponding to 
the Hamiltonian (fy) is given by 



1 d 2 ip(x) 

2 dx 2 



+ a(l - cos(x))ip(x) = Etp(x). (5) 



■0 is the wave function of the pendulum and E its energy. 
As the wave function of the pendulum has to be single- 
valued (0,0) , we impose the periodicity condition 



if)(x + 2tt) = tp(x). 



(6) 



It is possible to get an analytical solution of (g) by per- 
forming the following change of variables: u = (tt — x)/2 
and </>(«) = ip{x). We immediately obtain the canonical 
form of the Mathieu equation (|f47|) 



d 2 cf>{u) 



[a- 2qcos(2u))0(u) = 



where 



q = Aa and a = 8(E — a). 



(7) 



(8) 



Because of the previous change of variables, 4>(u) is now 
a 7r-periodic function. It can be shown ([|47|) that the 
Mathieu equation supports 7r-periodic solutions if and 
only if the characteristic value a belongs to an infinite 
countable set of values denoted by {a,2n{q), b2n(q)} n £ 
N. a,2n is related to the even Mathieu function ce2 n (u, q) 
whereas &2n is related to the odd Mathieu function 
se2n(u,q). As follows from their definition (||), q is di- 
rectly related to the energy of the separatrix E s = 2a 
that is to the depth of the cosine potential appearing in 
the Schrodinger equation whereas, up to a scaling and a 
shift factor, a represents the eigenenergies of the pendu- 
lum Hamiltonian. For a given value of q (that is of a), the 
series {a2 ra ,&2n} increases monotonically with n keeping 
the property &2« < »2n- For n = the only possible so- 
lution is an even function ceo(u,q) which represents the 
ground state of the system. 

According to the previous results, the analytical solu- 
tion of the Schrodinger equation (^) reads 

,AW - 1 re ( 7T ~ X n\ ■ - ° 2 "( g ) , „ 



for n&W and 



1 



-.ce Q 



-,q 



E 



(e) _ go(g) 



a (10) 



for n — which represents the ground state of the pen- 
dulum Hamiltonian. 



2. The Fourier representation (FR) 

According to the periodicity condition (||) imposed to 
the wave function, it is rather natural to consider its 
Fourier expansion: 



(ii) 
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where 



Vv, 



'2it 



ip(x)e 



c dx. 



(12) 



Because of the periodicity of the wave function, the 
Fourier space associated to the pendulum problem is in- 
finite but discrete. In the FR, the stationary Schrodinger 
equation becomes 



m 2 ip m - a(ip m+ i + ip m -i) = Eip v 



meZ (13) 



where for later convenience we have redefined the energy 
as 



E = 2{E-a). 



(14) 



Equation (13) represents a tight-binding equation 
whose hopping terms and onsite potential would be 
respectively —a and to 2 . In this equation, the sum 
("0 m +i + V> m -i) related to the cosine potential plays now 
the role of a discrete Laplacian while the onsite "po- 
tential" term m 2 comes from the kinetic energy. Thus, 
in FR, the terms of the Schrodinger equation have in- 
verted their role. This fact is important in understanding 
how the Discrete Wentzel-Kramers-Brillouin (DWKB) 
method ( 1 49 1 , f30f ) applies to the Mathieu equation (see 



section [n_B6|V 

Equation (|13|) is invariant under the transformation 
(to — > — m). This symmetry operation corresponds to the 
inversion of the momentum of the system (p — > —p) which 
originates in the time reversal symmetry of the original 
problem. This allows us to separate the eigenstates of 
( |l3| ) into symmetric |s) and antisymmetric \a) states. 

We may write the Hamiltonian of ( |l3| ) as 



H = m 2 |TO,)(m| — a(\m + 1)(to| + |to)(to + 1|) 



(15) 



where we have used the ket notation |to) to represent 
the plane wave function (x\m) = e imx /y/2w. The ma- 
trix representing H in the FR is infinite, tridiagonal and 
symmetric. Its diagonal elements are m 2 and the off- 
diagonal elements are constant and equal to —a. More- 
over, H m ^ n — iJ_ m ._„. This additional" central" symme- 
try is a direct consequence of the time reversal symmetry 
and the Hermitian properties of H . 



3. Low-energy states 

Without any onsite potential (a ~ 0), it follows from 
( |l5| ) that | to) and | — to) are eigenstates of H with iden- 
tical energy to 2 . Consequently, each state except the 
ground state |0) is two- fold degenerated in this limit. 
This two-fold degeneracy is due to the two equivalent 
possible motions consisting in rotating in a given sense 
or its opposite. 



The switching on of an onsite potential governed by the 
parameter a lifts this degeneracy but in a way depend- 
ing on the energy level of the state under consideration. 
Indeed, when a becomes nonzero, the pendulum system 
admits a new kind of motion, namely, the oscillating mo- 
tion which corresponds to the motion in the cosine arch 
of the potential well in the classical system (see Fig. [l]). 
The quantum system also admits such kind of states but 
their number is limited by the value of a because of the 
quantization rules. As a becomes larger, the number 
of states below the separatrix increases and can be esti- 
mated to 8yfa/ir by using the Weyl formula (see below 
( p7| ) in the limit k — > 1). If the value of a is large enough 
it becomes possible to look at the properties of the low 
energy states by expanding the cosine potential around 
zero in a series of powers of x. We obtain 



a(l — cos(ie)) = a 



(16) 



The first order of this expansion leads to the har- 
monic approximation of (Q) whose frequency is given by 
lu = s/a. Thus, low energy states of the quantum pendu- 
lum are well represented by the corresponding low energy 
harmonic eigenstates (at least for large enough values of 
a). As the level spacing of the harmonic oscillator is 
constant (and equal to u>) we expect the states to be 
quite regularly spaced deep inside the well. This is il- 
lustrated in Fig. [I] which represents the lower part of 
the spectrum of a pendulum whose parameter a = 50. 
The energy levels of \s) and \a) states have been respec- 
tively represented to the right and to the left. We clearly 
see that the few first states are quite regularly spaced 
while they condensate in reaching energies close to the 
separatrix. By using standard second order perturbation 
theory within the natural harmonic basis, we may obtain 
the first corrections to the pure harmonic spectrum. This 
perturbation in the nonlinear terms gives 



En = ^(n+(n+l))u-±-(n 2 + (n + l) 2 ) 
- ^(" 3 + (- + l) 3 ) 1 + ofA 



(17) 



where n labels the states (n £ K). Of course, this expres- 
sion ceases to be valid when nonlinear corrections become 
large as compared to the harmonic term. Nevertheless, 
it shows how the spacing becomes smaller as the energy 
increases but remains below the energy of the separatrix: 
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FIG. 1: Energy levels of a pendulum with a = 50. States 
have been distinguished according to their parity. Only a few 
states lying above the separatrix have been displayed. The 
number n of the corresponding doublet is written on top of 
it. 



4- High-energy states 

Above the separatrix, Fig. [I] shows that symmetric 
and antisymmetric states glue together to form pairs of 
nearly degenerated states. From a physical point of view 
this phenomenon is due to the fact that far above from 
the separatrix the cosine potential appears like a pertur- 
bation of the free rotor and only shifts the levels a bit 



around their value £: f i; corotor 



m 



(to 3> 2a). Because 



of the momentum reversal symmetry, the Hamiltonian 
( |l5| ) diagonalizes into two blocks each related respectively 
to \s) or | a) states. In the FR, the corresponding reduced 
matrices are still tridiagonal. By applying standard per- 
turbation theory in the parameter a to one of these ma- 
trices whose spectrum is now free of nearly degenerated 
eigenvalues, we obtain corrections to the free rotor energy 



E 



free rotor 



m 



E, 



To second order, we find 



2a 2 



An 2 - 1 



O 



(19) 



Explicit corrections up to a 6 /n 10 are obtainable (see |47j 
p. 724). As we mentioned, this energy is the mean energy 
of a doublet. This is indicated by the subscript d in (19). 



blet. Because of their symmetry, the expectation value 
of the momentum for \s) and \a) are zero. 

A possible way to get this splitting is to use high order 
perturbation theory as it has already been done in |l2]| 
(equation (12)). A proof of this formula has been derived 
in Q and applied to the quantum discrete self-trapping 
equation. The special " centro-symmetric" form of the 
tridiagonal Hamiltonian matrix in the FR allows us to 
compute the exact leading order of the splitting as 



AE n = ~A£„ = Ua 2n {Aa) - 6 2n (4a)) 
A o 



,.2n 



(2n- 



l)! 2 ' 
(20) 



Notice that due to the obvious relation between the eigen- 
functions of the SPP and the Mathieu functions, the 
splitting AE n gives in turn the result of the splitting 
between the characteristic values associated to the sym- 
metric and antisymmetric 7r-periodic Mathieu functions 
ce2n and sei n . A similar computation allows us to obtain 
the splitting for the 27r-periodic Mathieu functions. By 
using the notations of p7[ , we get 



a r (q) - b r (q) 



(r-l)\ 2 



(I)' 



o Iq 



r+n 



as 



(21) 



To our knowledge this is the first time that such a split- 
ting has been exactly computed in leading order p8| . In 
particular, it corrects and precises the formula obtained 
in p7| . Moreover, it can be verified to be correct by us- 
ing the small q expansions of the low order characteristic 
values themselves. 

As n, which represents the label of the doublet, be- 
comes large enough, we find the following asymptotic 
form for the splitting 



AE r , 



1 

127T 



2/7 



(n 



(22) 



indicating that it decays more than exponentially fast as 
n increases. Moreover this expression contains its limit 
of validity as the number n has to be greater than the 
critical value n c — e-y/q/2 to ensure that the splitting is 
small (see comment ]69[). 



6. Discrete WKB theory 



5. Analytical computation of the splitting: High-order 
perturbation theory 

It is also possible to calculate explicitely the splitting 
between the symmetric and antisymmetric states of a 
doublet labeled by to (to 2 ^S> 2a). The inverse of this 
splitting is a direct measure of the time required by a ro- 
tating state to invert its initial momentum. It has to be 
noticed that by "rotating state" we mean the superposi- 
tion of the states |s) and |a) belonging to the same dou- 



Another way to compute the doublet splitting is to 
use the DWKB method developed in [[l9) and to use the 
discrete counter par t of the Herring's formula as done in 
f30f . Following |5(J and according to (|l3|), the splitting 
between the nearly degenerated states reads 

AE n = \\E^-E^\= a s^a^ (23) 

where and a± are the Fourier components of the 
symmetric and antisymmetric states of the n-th doublet. 
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Instead of evaluating these components by using the con- 
nection formulae of j49|, we use directly the resulting 
forms given in |47j]. The asymptotic behavior (see ex- 
pression 20.2.29 of ^7j) of the Fourier components of 
the Mathicu functions yields 



4 2n) (9) 



[n 



l)! 



q" 



n!(2n- 1)! 4 



2n(<?) 



A 2n 
<1 



(n- l)!(2n- 1)!4«-! 2nW 



where A and B are the coefficients of the cosine and sine 
Fourier series of the 7r-periodic Mathieu functions cem 
and se2n- The superscript indicates the order whereas 
the subscript labels the Fourier components. The cor- 
rect normalization yields Sg™' = \f2A^ n ^ (q) and = 

B^ n \q)/\/2 (see @ and (|l|)). By taking into account 
that A 2 2 l{q) = B^(q) = 1 + o{q) and a = q/A, we finally 
get 



,.2n 



[(2n- l)!] 2 



(24) 



which coincides with (|20|). The graph (Fig. g) shows 
the spacings and the splittings of the SPP spectrum for 
a = 50 as a function of the energy E. The crosses cor- 
respond to E n+ \ — E n as obtained from the numerical 
diagonalization of (|l3|). The two upper branches (dotted 
line) represent the spacing between neighboring eigenval- 
ues (when E < E s ) and the spacing between the doublets 
(when E > E s ) computed by using the density of states 
(p6|). They are shown to be in excellent agreement with 
the numerical data. The splitting, that is, the energy 
difference between the states \s) and \a) of the doublets 
is represented by the third branch which is decreasing 
rapidly as the energy increases. Again, the analytical 
result (solid curve) given by ([24]) is excellent. The in- 
set shows a comparison between this analytical splitting 
and the numerical result of the numerical splitting ob- 
tained by using Herring's formula (E3rj. Notice that the 
splitting can be now as small as lO^ 1 ^ although it has 
been computed with a simple FORTRAN scheme in dou- 
ble precision (16 digits). Because the eigenvalues are of 
order of unity it is of course impossible to obtain such 
a result by subtracting two neighboring eigenvalues ob- 
tained by diagonalizing ( |l3| ) in double precision. The 
precision is in this case limited to ~ 10 -14 . But the way 
to compute the eigenvectors makes it possible inasmuch 
as the numerical limit becomes not the number of digits 
but the smallest number representable by the computer. 
This result shows that it is possible to compute very small 
splittings by using common FORTRAN routines instead 
of using high precision schemes provided by Mathematica 
or MAPLE for instance. 

An equally important property of the DWKB method 
is that it provides with a simple and nice picture of the 
basic properties of the eigenstates. Indeed, as explained 
in p9| , it is possible to associate a "classical" Hamilto- 
nian defined by H c i = m 2 — 2a cos 4> to the three-term 





.■It H — r— 




4 1 


Numerical Splitting 
Analytical Mean Spacing 
Analytical Pair Splitting — 


f 


\ 1B-20 

4 ra 1e " 40 
\ I le-60 
\ % 1B-80 

T 1B-100 

\ 16-120 


%L Herring + 

- ^Su Analytical - 




200 400 600 800 1000 1200 
Doublet energy 



300 

Energy 



FIG. 2: Spacings and splittings of the spectrum of a quantum 
pendulum for a = 50. Quantities are described in the text. 



recursion relation ( |l3| ) , where m and <j> correspond to the 
conjugated "coordinate" and "momentum". This defi- 
nition is rather natural as <f> = —id/dm is precisely the 
Fourier representation of the angular variable x. Inter- 
preting as a momentum, the expression of H c i shows 
that its classical motion is confined between two "po- 
tential" curves defined by U (m) = m 2 ± 2a (see Fig. 
@). The idea is then to use this fact to compute the 
DWKB solution related to this "classical" Hamiltonian. 
The internal region defined by the two parabolas U (m) 
represents the classical allowed region whereas the re- 
gions outside are forbidden. The eigenfunctions are thus 
localized on the allowed region while they decay expo- 
nentially (or faster) in the forbidden ones. This provides 
a natural way to express their location and their exten- 
sion on the Fourier lattice according to their energy. We 
find indeed that at sufficiently high energy (far above the 
separatrix) the eigcnfunction of energy E ~ m 2 /2 + a is 
localized around m and — m whatever its parity. Its local- 
ization length around these two centers is roughly given 

by L(E) = \fW- 



2a - \f~E- 2a - E s /y/2(E-a) as 
E 3> E s = 2a. The localization length in the Fourier 
space thus decreases like the square root of the energy. 



7. Density and number of states 



We end this section devoted to the properties of the 
single pendulum problem by giving the expressions of the 
number n(E) and the density p(E) of states computed 
by the mean of Weyl formula: 



P(E) 



1 

27 



5{E - H(x,p))dxdp. 



(25) 
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m : Fourier components 

FIG. 3: "Potential" curves U ± used in the DWKB theory 
to compute the eigenfunctions of the SPP. Some symmetric 
(solid) and antisymmetric (dashed) eigenfunctions have been 
superposed to them. Their base line situated at their energy 
level allows to compare the location and the width of their 
peaks to the interval between the potential curves U ± . The 
value of a is 50. For the sake of visibility, the wave functions 
have been multiplied by a suitable scaling factor. 



It yields 



P(E) 



^LktC(k) , E < 2a ; k — \[^- 
ttVE V 2a 



2^2 



ttVE 
and by integrating 



E > 2a ; k= \ — (26) 
E 



n(E) = 



8^a~ 



(£(k) - k' 2 K{k)) , E < 2a 
£{k) , E > 2a (27) 



with the same definition of k than above. In these ex- 
pressions £ (k) and K{k) denote the usual complete ellip- 
tic integrals of the first and second kind (see e.g. pTf ). 
The above expressions of p(E) show that the density 
of states develops a logarithmic divergency close to the 
separatrix. In this limit, its expression reads p(E) ~ 
\n(16E/\E — E s \)/ (ic-s/a). This phenomenon is known as 
a Van Hove singularity J5l| and has also been discussed 
for the case of the quantum dimer JT^ ] . The graph (Fig. 
|J) shows the excellent agreement between the density of 
states (DOS) computed by the mean of Weyl's formula 
and the DOS computed numerically by diagonalization 
of©. 

Using (|27|), we may now obtain the splitting AE as a 
function of the doublet energy E: 



AE(E) 



(E 



(28) 
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FIG. 4: Density of states for a pendulum with parameter 
a = 5. 10 3 . Numerical data have been obtained from ( Jl5| ) for 
a matrix truncated to 1001 Fourier components (—500 < m < 
500). Comparison with analytical expression (E6|) shows an 
excellent agreement until truncation errors become important 
(E ~ 1.25 10 s ). 



III. TWO COUPLED PENDULA 



The Hamiltonian of the two coupled pendula is given 



by 

H = H!+H 2 +Hi 



lilt 



2 

E 

1=1 



+ a(l — cosxi) 



+ e(l - cos(xi - x 2 )) 
(29) 



where e > represents the coupling parameter between 
the pendula. 

The coupling has been chosen to be periodic in {x\— x%) 
to allow solutions where one pendulum is oscillating 
whereas the other one is rotating. This condition is of 
course essential to obtain rotobreather type solutions. 
Notice finally that such a Hamiltonian has already been 
used to describe classical nonlinear rotating modes in a 
chain of coupled pendula J52[. 

The Hamiltonian ( ]29j ) possesses two different symme- 
tries. One is related to the exchange of the coordinates of 
the two pendula {S p : {x\,p\) «-* {x2,P2)} and the other 
to the reversal of the global momentum of the system, 
{S m : (jpi + P2) «-> —{pi+P2)}- The latter represents 
the generalization of the momentum reversal symmetry 
already observed in the single pendulum system. 



A. Classical rotobreather 

Although the system is invariant with respect to the 
permutation of the coordinates, some solutions of the 
equations of motion are not. Indeed, there exist exact 
solutions of the Hamilton equations derived from (|2£ 
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which consist of two pendula oscillating at the same fre- 
quency but with different amplitudes. It is also possi- 
ble to obtain solutions where one pendulum is oscillating 
whereas the other one is rotating, both again evolving at 
the same frequency. Hereafter, both kinds of solutions 
will be respectively referred to as "breather" and "roto- 
breather" solutions. Of course, because of the size of our 
system, the exponential spatial decay property of usual 
breathers becomes meaningless and the classical breather 
type solution refers only to exact time-periodic solutions 
which break the permutation symmetry S p . 
As Poincare sections (PS) are known to be a useful tool in 
describing the behavior of non integrable dynamical sys- 
tems, we will use them to locate the orbits of the classical 
rotobreathers at a given energy. This energy has to be 
larger than the separatrix level of the SPP to allow one of 
the pendula to rotate whereas the other one is at rest in 
the uncoupled (anti-continuous) limit (e = 0). But this 
energy can not be too close to the SPP separatrix as e 
becomes non zero. Indeed, the system develops a chaotic 
layer in its vicinity and prevents any stable periodic solu- 
tion of the rotobreather type from existing. Concerning 
the classical rotobreather itself it is obtained by using 
standard numerical methods like a Newton scheme or 
a variational method (see e.g. p5|,p9]]). The purpose 
of this paper is not to carry out an extensive study of 
the classical system but merely to deal with its quan- 
tum counterpart. Consequently we will just give here an 
example of a classical rotobreather. The values of the 
parameters are a = 5, e = 1. The energy of the roto- 
breather is E b ~ 31.77 and its period T b ~ 0.885. The 
initial data are p\ ~ 7.969, p\ ~ -0.160, x\ = and 
x\ = —5.2 10~ 9 . The first pendulum is rotating whereas 
the second one is librating. The graph Fig.|| represents 
two Poincare sections realized with the condition X\ = 
at an energy fixed to E b . The first PS is plotted in the 
phase space of the second pendulum. As its trajectory is 
strictly periodic, the Poincare section of the rotobreather 
orbit simply consists of a single point. It has been rep- 
resented by a "star" symbol. Other periodic trajectories 
present in this system have been represented by a " cross" 
and a "times". They respectively show the in-phase and 
out-of-phase motions which are symmetric modes and 
thus can not be considered as breathers. The second 
PS represents the same situation in the momenta space. 
Given the natural Fourier representation of the quantum 
problem due to the 27r-periodicity of the wave function, 
it provides the ideal frame for the comparison of the clas- 
sical and the quantum situations. Notice that because of 
the existence of symmetries the presented rotobreather 
solution is not the only one. Indeed, in the momentum 
space (pi,p2), the global momentum reversal symmetry 
is nothing but the reflection symmetry with respect to 
the point (pi — 0,p2 = 0). And the permutation sym- 
metry is the mirror symmetry with respect to the line 
Pi = P2- Using these two symmetries we easily obtain 
the location of the four classical rotobreather solutions 
(±Pi,±p!)- Nevertheless, and despite these symmetries, 



if the system starts initially on one of the classical ro- 
tobreather orbits, it remains on this orbit for an infinite 
time. This prevents the classical system from transferring 
the initial excitation from one pendulum to the other. We 
will see that such a situation is impossible to be realized 
in the quantum system. 

In order to find the quantum counterpart of the clas- 
sical rotobreather, we will proceed exactly as in the clas- 
sical case. Starting from the anti-continuous limit where 
we know the classical rotobreather to exist, we will find 
the corresponding quantum state and study its evolution 
with respect to the coupling e. 

B. 2D Fourier space 

The state of the two pendula system tJ){x\,X2) has to 
be 27r-periodic in each of its variables as in the single pen- 
dulum problem. We thus expand it as a double Fourier 
series 

\ip) = X! ^ m ,n\m,n) (30) 

(m,n)GZ 2 

where Wn,n) = \m) <8> \n) and (x\n) — exp(inx)/y2n. 
With (p9|), the eigenvalue equation H\ip) = E\tp) reads 

E-ipm.n = {m 2 + n 2 )i> m . n - a(V>m+l,n + i>m-l,n) 
—a(ip l) - £(-i/Wl,n-l + Vm-l.n+l) 

(31) 

where E = 2(E — 2a — e) is a shifted and rescaled energy. 
In this discrete Fourier space, the symmetry operations 
S m and S p become: {S m \{m,n) — > (— to, — n)} (mo- 
mentum reversal) and {S p \(m,n) — > (n, m)} (permuta- 
tion) . As H is invariant under these symmetries, it com- 
mutes with the corresponding operators and its eigen- 
states gather naturally in four different symmetry classes. 
We denote the four different kinds of eigenstates by 
{\s}, | a), \s), | a)}. Their symmetry properties are listed 
hereafter. 



State 


Momentum reversal S m 


Permutation S p 


\s) 


S 


^m,n — 3 — m, — n 


S 


°m,n — $n,m 


\a) 


A 


®m,n — Q—m. — n 


s 


@"m,n — Q"n.m 


\s) 


S 


Sm,n = S— m, — n 


A 


Sm,7i — Sn.m 


\a) 


A 


^m,n — & — m,— n 


A 





S and A indicate respectively that the state is symmet- 
ric or antisymmetric with respect to the corresponding 
symmetry. To evaluate numerically the eigenvalues and 
eigenvectors of ([^), we truncate the (normally infinite) 
system to — N < m,n < N, where N is chosen suffi- 
ciently large to prevent any important truncation errors 
for states whose energy E <C 2N 2 . In this case, the total 
number of computed eigenvalues is (2N+ 1) 2 . Due to the 
symmetries, we only diagonalize the sub-matrices repre- 
senting the four different classes of eigenvectors. The 
rank of the sub-matrices related to |s), \a), \a), and \s) is 
(N + l) 2 , N{N + 1), N(N + 1) and N 2 respectively. 
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FIG. 5: Poincare sections of the two coupled pendula system in the phase space of the second pendulum (left) and in the 
momenta space (right). The conditions of the maps are X\ — 0, pi > (left) and xi = 0, X2 > (right). The energy E = 31.77. 
In-phase (cross), out-of-phase (times) and rotobreather (star) trajectories are represented. 



The global spectrum at the anti-continuous 
limit (e = 0) 



The corresponding eigenstates are given by a product of 
correctly symmetrized plane waves. Namely, 



There are two limiting cases where the eigenvalue equa- 
tion ( |3l| ) can be solved analytically. These correspond to 
situations where the classical system becomes integrable. 
The first one is realized when the coupling parameter e 
is equal to zero, i.e. where the system consists of two 
identical decoupled pendula. The second one is realized 
when the onsite parameter a becomes zero, where, by 
passing to the center-of-mass representation, the system 
can be reduced to a free rotor plus a decoupled pendulum. 
In these two limits the global spectrum of the system is 
given by the sum of two one-particle spectra. Moreover, 
in the limit where e is equal to zero, the two spectra are 
identical. This leads to a two-fold degeneracy of the main 
part of the global spectrum. Nevertheless, by construct- 
ing the eigenvectors in such a way that they belong to 
a given symmetry class, the eigenstates are unambigu- 
ously defined and already represent the proper zeroth 
order states suitable for any perturbation calculation in 
e. 





1 








s) 


= 2 
1 


\m 


n) + 


— to, 


a) 


= 2 


\m 


n) — 


- m, 




1 








s) 


= 2 


\m 


n) + 


- m, 




1 








a) 


= 2 


\m 


n) — 


— m, 



\n, m) 
In, m) 



n, —to)) 
n, —to)) . 



Notice that E m n is a nonnegative integer in this limit. 
Its degeneracy can be expressed as follows (result due to 
Gauss, see e.g. 53 1): 
Let E e N* and 

g[E) = #{(to, n) E Z 2 : m 2 + n 2 = E} 
be the degeneracy of the level E. Let 



n 



E = | | p 

prime p 



e(p) 



1. Limit of zero onsite parameter: (a = 0) 



be its prime number factorization. If V p = 3 mod 4, 
e(p) = mod 2 , then E can be written as the sum 
of two squares, i.e. g(E) ^ (Fermat). In this case 
(Gauss), its degeneracy is given by 



Let us first derive some results concerning the anti- 
continuous limit (e = 0) when the onsite parameter is 
itself equal to zero (a = 0). In this case, the system 
consists of two free rotors and the global spectrum is 
straightforwardly given by 



E m .n = m 2 + n 2 . 



(32) 



9(E) = A J] «P) + 1 )- 

p=l mod 4 



(33) 



Moreover, the number of states (i.e. the integrated den- 
sity of states) behaves asymptotically as 



i(E) = J2 9{E') = nE + 



(£-►00). (34) 



E'<E 
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2. Switching on the onsite potential (a 7^ 0) 

The last expression is not only valid in the limit of 
free rotors but also when a and e are non zero. This 
is confirmed by the Weyl's formula which indicates the 
zeroth order term of the degeneracy g(E). Indeed (Q 
p.497, 21.4) 

gwcyi(E) = - ( — / dxi dx 2 I (35) 

* \^ Jv(x 1 ,x 2 )<E J 

where —V(xi, x-i)j2 = a(cosxi + COSX2) +s cos(a;i — x^)- 
The prefactor 1/2 is coming from the rescaling of E as 
compared to E. As (xi, xz) s] — n, ir] 2 , when E > E + = 
max (V(xi, X2)), the integral term is constant and equal 
to (27r) 2 . Thus gwcyi(E) — ir and nw ey \(E) = ttE + const 
when E > E + . It is possible to show by direct calculation 
that the constant term is equal to zero and that E + = 
2(2a - e) if a > 2e or E+ = 2(e + a 2 /2e) if a < 2e. As 
E + gives the " top" of the potential for any values of the 
parameters e and a, it also represents a natural transition 
value in the global spectrum indicating the level at which 
the influence of the potential V starts to be weak. 

Expression ([mJ) thus provides a very useful guide for 
checking the energy level at which truncation errors be- 
come important. Note that truncation errors are of two 
types. First, because of the truncation of the whole ma- 
trix, there exists an energy threshold for which some 
eigenvalues of the spectrum are missing. Indeed for e = 0, 
the global spectrum is given by the sum of the SPP spec- 
trum with itself. In this case, even assuming that the nu- 
merically computed eigenvalues are exact, the SPP spec- 
trum ends at E° pp because of truncation. Let us denote 
the energy of its ground state by Eq P . Then the first 
eigenvalue to be missed is Ef p f + Eq PP where Ef p f is 
the first eigenvalue following i?f pp . This provides with 
a threshold where the eigenvalues start to be ranked in 
a wrong way. Their label becomes false and so does the 
computed number of states n(E). This is shown in Fig. 
^ which represents the spectrum of two uncoupled pen- 
dula as computed from equation ( |31~j ) for different values 
of N (recall that — N < m,n < N). The second type of 
error induced by the truncation is a modification of the 
values of the energies themselves. As we have verified 
these errors increase as we reach the upper end of the 
truncated spectrum but are nevertheless very small if we 
respect the threshold indicated above. 

Although the degeneracy (|3^) has been obtained for 
the case where both the onsite and the coupling param- 
eters are zero, it provides with useful information when 
considering the near degeneracy of n-uplets arising at suf- 
ficiently high energy (that is, far above twice the energy 
of the separatrix of the SPP). 

As the onsite parameter a becomes different from zero, 
the first part of the free rotor spectrum of the SPP is 
modified. Far below the SPP's separatrix, the quadratic 
spectrum Ef FP =l 2 is replaced by a harmonic oscillator 
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FIG. 6: Spectrum of the global uncoupled (e = 0) system 
computed by diagonalizing ([H]) for different values of the 
maximal number of Fourier components iV of the wave func- 
tion The linear dimension of the Hamiltonian matrix is 
(2iV + l) 2 . q = 50. E + (see text) has been indicated as 
well as the thresholds where, because of truncation, the com- 
puted spectrum starts to be larger gapped. Thresholds are 
represented by straight lines ended by the symbol of the cor- 
reponding spectrum. The asymptote E — n/ir is also shown. 



one i?f pp ~ 2(i/aZ — a). The number of states n(E) of 
the global uncoupled system thus becomes a quadratic 
function of the energy. This explains the form of the 
function E(n) observed in Fig|| which starts as a square 
root and asymptotically becomes a linear function of n. 



D. Quantum rotobreather state 

1. Anti- continuous limit 

The goal of this paper is to define and to study a quan- 
tum rotobreather state |$) whose properties are very sim- 
ilar to those of the classical rotobreather. Consequently 
it is natural to look for a quantum state which, at the 
anti-continuous limit [e = 0), represents a state consist- 
ing of a rotating pendulum and another one at rest (up to 
quantum fluctuations). As the pendula are not coupled, 
this state is represented by the tensorial product of the 
SPP's ground state |ito) (corresponding to the pendulum 
at rest) and a superposition of two states (|c n ), \a n }) be- 
longing to the same doublet n of the SPP (correspond- 
ing to a rotating pendulum). The addition of the two 
SPP spectra and the construction of the quantum roto- 
breather |$) is schematically depicted in Fig. [?]. The 
states denoted by \a) (\a)) are symmetric (antisymmet- 
ric) with respect to the inversion of the momentum: 

1$) = 4=ko) ® (|<Tn> + K»- (36) 
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FIG. 7: Schematic representation of the sum of two pendula 
spectra. Straight solid arrows indicate the levels to be added 
and dashed arrows the symmetric (permutation) operation. 
The result is indicated in the global spectrum by a curved 
arrow. The construction of the quantum rotobreather state 
|$) is explicitely represented. 



In order to obtain the state |<£>) at e = as a linear 
combination of the corre ctly symmetrized eigenvectors 
listed in the section 1IIB, we first construct it in terms 




of the SPP eigenvectors. The correct symmetrization of 
the tensorial products leads to the following expression 
of the four possible classes of states: 



or 



or 



T2 iWU 
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+ \ct) 
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-\oi) 
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5|a) 
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5|a) 


- \a) d 


i\a)) 



(37) 



where \a), \a J ) (|a), |a')) represent any of the symmetric 
(antisymmetric) SPP eigenvectors. 

Using these expressions, |$) is readily written as 



1$) 



(l*> 



(38) 



where \s),\s) = [\a ) ® \<r n ) ± Wn) ® \<r ))/V2 and 
|o),|o) = (|cr ) ® |a„) ± |ct„) <g> |a ))/\/2. In this anti- 
continuous limit, we may compute the eigenvalues corre- 
sponding to each of the four states which contribute to 



|$). We obtain 

E s = Eg = E ao + E„ n ; E a = E a = E aa + E an . (39) 

As the expression ( p^ ) shows, the absence of any coupling 
between the pendula is responsible for the true degener- 
acy of the states whose parity is the same with respect 
to the momentum reversal symmetry . This follows from 
the fact that the global spectrum is given by the sum 
of the SPP spectrum with itself in the anti-continuous 
limit (cf. Fig. 0). Thus the energy level E ao + E an is 
obtained in this limit either by adding the energy of the 
ground state \<jo) of the first pendulum to the energy of 
the symmetric state \a n ) of the n th doublet of the second 
pendulum, or by adding the energy of the ground state 
of the second pendulum to the energy of the symmetric 
state of the first one. 

We have already shown (see |20|) that the presence of an 
onsite potential (a ^ 0) lifts the degeneracy of the sym- 
metric and antisymmetric states belonging to the same 
doublet of the SPP. We may thus compute the splitting 
between the two pairs of degenerated states (|s), \ s)) and 
(|o),|o»: 



,.2n 



E R — E„ 



(2n- l)! 2 



(40) 



For sufficiently large n this splitting is extremely small 
and thus the four states making up |$) form a quadru- 
plet of nearly degenerated eigenenergies which has been 
represented in Fig. 0. Notice finally that the inverse 
of this splitting is directly related to the time taken by 
the system whose initial state is |<E>) to reverse its total 
momentum. As one of the pendula is at rest (up to quan- 
tum fluctuations), this means that the pendulum initially 
rotating in a given sense tunnels into the state which cor- 
responds to the opposite rotation sense on a time scale 
t p ~ l/5 n . This effect is purely quantum. 



2. Nonzero coupling e 

As soon as the coupling parameter e is nonzero the 
degeneracy of the states (|s), \s)) and (|a), \a)) of the 
quadruplet is lifted. Here we are interested in the com- 
putation of the corresponding eigenvalues in order to find 
the splittings which determine the evolution of quantities 
such as individual energies or momenta of the pendula. If 
the splittings between the states making up the quantum 
breather become of the order of the mean level spacing 
of the spectrum, we may conclude that the breather so- 
lution is lost |5q| . 

One possible source of a dramatic increase of the split- 
tings can be a strong overlap of the quantum state with 
the chaotic layer of the classical system (see e.g. Q for 
a general review of the manifestations of classical phase 
space structures in quantum mechanics) . This overlap is 
generally computed by the mean of a Husimi distribution 
which is one of the possible phase-space representations 
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of a quantum state (see e.g. [157]]). This distribution is 
then superposed to the corresponding Poincare section 
of the classical system which allows to compute the over- 
lap. This method has for instance been used in the case 
of a driven bistable system |5j| . At the same time these 
studies have shown that doublet states overlapping up to 
70% with the chaotic layer may still possess a small split- 
ting. To avoid misinterpretations, here we will compare 
the Fourier components of the quantum rotobreather |$), 
with the Poincare section of the classical system in the 
momenta space. The main result will be that the phase 
space locations of the classical and the quantum roto- 
breathers are roughly the same. 

This direct comparison is possible due to the results 
obtained in |59 showing that the Husimi distribution of 
the pendulum problem can be found analytically. It fol- 
lows that for the pendulum potential, the discrete Fourier 
representation of the eigenfunctions and their Husimi dis- 
tributions restricted to the momentum space differ in- 
significantly. 



3. Tunneling of the |<3?) state 

In this section we derive the expressions of some rel- 
evant quantities which allow us to follow the time evo- 
lution of the initial state |$) = (|s) + \s) + \a) + |a))/2. 
This state is formed at e ^ by the eigenstates which be- 
long to the quadruplet of states. These states in turn are 
defined by the tensorial product of the ground-state of 
the SPP spectrum and one of its doublets at the anti- 
continuous limit. Because |<I>) is not an eigenstate it 
evolves in time. In order to visualize this evolution and 
because we are working in the 2D discrete Fourier space 
(2DFS), we may compute the time evolution of its mo- 
mentum. This is done by computing the two functions 
(AP) = ($|A - P 2 |$) and (P) = ($|A + P 2 |*) where 
Pi and P2 are respectively the momenta operators of the 
pendula 1 and 2. If the total momentum of the system 
were conserved (as in the integrable limit a = 0), the 
difference of the momenta (AP) would exactly represent 
the transfer of momentum between the pendula. This is 
not the case for a nonzero value of a. Although the total 
momentum (P) is in general not conserved, for not too 
small values of e it evolves very slowly as compared to 
(AP). We thus may consider that the latter represents 
the transfer of the excitation from pendulum to pendu- 
lum with a good accuracy. 

Another quantity of interest is the difference between the 
individual (or onsite) energies (AP) = ($|Pi — P2I*) 
where Hi = Pf /2 + a(l — cos JSQ). As the total energy of 
the system is conserved, (AP) measures the transfer of 
energy between the pendula. 

In order to give a simple expression of the above quan- 
tities let us use the two symmetries of permutation and 
momentum reversal. We denote their respective opera- 
tors by V and M.. It is easily shown that they are uni- 



following relations hold: 

M^PiM = -P 
M ] HiM = iP 
M^H mt M = Pint 



V ] P{P = P 



P T P ? P = H 



j 



ptPintP = Pi, 



(41) 
(42) 
(43) 



where G (l,2),i =/= j and where Pint 

cos(A"i — X2)) is the interaction energy term. 
Using the fact that the Hilbert space £ associated to 
the coupled pendula problem can be written as the di- 
rect sum £ = £ s (B £g ® £ a ® £ a , any operator may 
be represented by a 4 x 4 block-matrix in the basis 
({l s i)}j {|%)}j {l a fc)}i {l^z)}) formed by all the states be- 
longing to £ s first, then all the states belonging to £ s , 
etc... in the order described above. From the preceeding 
relations and by using the symmetry properties of the 
eigenstates, it follows that the relevant operators for our 
study have the form 



P. 



( 





P 

1 as 


±Pas 








±Pas 


P— 
J as 


p 

1 sa 


±Psa 








\±P.a 


P— 
J sa 









(44) 



and 



Hi, 2 = 



( H ss 
±H S 


V 



±Pss 



















H aa 


±Paa 





±Paa 


H aa 



(45) 



We have used the condensed notation M „ = (/i|Oi|f) 
which symbolically represents all the matrix elements of 
the observable 0\ between states of two given subspaces 
£ jl and £ v . The sign ± is + for 1 and — for 2. In this 
basis, the operator P; n t is block-diagonal. 
Using the time evolution of the initial state |$) 



\m) = \Y, e ~ iEi ' 



we finally obtain 



M *|/z) where (i<E{s,s,a, a} (46) 



(AP) = 


<s|^i 


a) 


cos(P s — 


E- a )t 


+ 


(s\Pi 


a) 


cos(P g — 


E a )t 


(P) = 


(s\Pi 


a) 


cos(P s — 


E a )t 


+ 


(s\Pi 


a) 


cos(Pj — 


E- a )t 


(AP) = 


(s\Hi 


\s) 


cos(P s — 


E s )t 


+ 


{a\Hi 


\a 


) cos(P a - 


-E- a )t 



tary: V^V = M^M = 1 (see e.g. |60|). Moreover, the 



(47) 
(48) 
(49) 



Notice that the energy differences occurring in (AP) are 
zero at the anti-continuous limit. No transfer of energy 
between the pendula occurs due to the fact that they 
are decoupled. The excitation (and thus the energy) is 
entirely conserved on its initial site. In this limit we 
know (in leading order) the value of the splittings E s — 
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E a = Eg — E a — a 2n /(2n — l)! 2 and thus the time taken 
by the system to reverse its initial momentum: r p = 
7r(2n— l)! 2 /a 2n where n labels the doublet of the rotating 
pendulum. These splittings also occur in (AP). Thus 
for small e we can not assign the meaning of a transfer 
of momentum between the pendula to the time evolution 
of(AP). 

In the other integrable limit where e ^ but a = 0, 
it is also possible to compute the exact spectrum of the 
system. Notice that the total momentum of the system 
is strictly conserved in this case. The absence of the 
onsite potential is also responsible for the degeneracy of 
the states of different parity concerning the momentum 
reversal symmetry. We thus have E s = E a and Eg = E a - 
It is possible to show (see appendix A) that the splitting 
E s — E s = (a n (—2e) — b n (— 2e))/4 which in leading order 
gives \E S - Eg\ = 2(e/2) n /(n - l)! 2 where n is the label 
of the doublet of the initially rotating pendulum. The 
time taken by the system to transfer its excitation from 
one site to the second one is r e = ir(n — l)! 2 /2(e/2)™ in 
this limit. By comparing r e with r p we observe that for 
large enough n, r p ^> r e holds, regardless the values of e 
and a. This indicates that in general transfer of energy 
(and momentum) from site to site will be a much faster 
process than the process of total momentum reversal. 

When a and e are both nonzero, the different split- 
tings have to be calculated numerically by following the 
evolution of the quadruplet as a function of the coupling 
e. This is the purpose of the following section. 



4- Following the quadruplet 




FIG. 8: Evolution of a part of the spectrum as a function of 
the coupling e. The onsite parameter is a = 5. Energies are 
rescaled according to E — 2{E — 2a — e). The quadruplet 
under investigation starts at an energy E ~ 41.43. It corre- 
sponds to a state formed of the ground-state of the SPP and 
the doublet n = 7. 

In Fig. ||, we have plotted the evolution of a part of 
the global spectrum of the coupled pendula as a function 



of the coupling e. The onsite parameter is a — 5. The 
energies of the different classes of states have been rep- 
resented by different lines and symbols. The quadruplet 
which represents the levels of the states {|s), \s), \a), \a}} 
from which |<&) is formed, has an energy E ~ 41.43 at 
e = 0. This corresponds to a quadruplet made of the 
ground-state of the SPP and the doublet n = 7 at the 
anti-continuous limit. Indeed, by using ( ^7|) and (|is|), 
we obtain E% ~ 49 + 50/195 and E ~ i/(5) — 1/16 => 
E Q ~ y/{5) - 1/16-10. The sum of these energies gives 
E ~ 41.43. At e = 1 the quadruplet energy is E ~ 41.54. 
This corresponds to E ~ 31.77 which precisely represents 
the ene rgy o f the classical rotobreather presented in the 
section III A. 



At slightly larger energy in the spectrum we observe 
an octuplet. It consists of states mixing two neighboring 
doublets, namely n = 4 and n' = 5. Its energy is around 
E ~ 42.3 which can be computed also by using (|lS|). 
The reason for the occurence of this octuplet is depicted 
in Fig. 0. By combining the two states of the doublet n 
with those of n' we obtain four different (although nearly 
degenerated) levels of energy. But there are two ways to 
attribute these doublets to the pendula 1 and 2 (repre- 
sented by the solid and dashed arrows in Fig. 0). We 
thus obtain two identical quadruplets which yield an oc- 
tuplet. Finally, the visible cluster of levels at E w 43.25 
corresponds to a quadruplet made of the combination of 
a single symmetric state (the last one situated below the 
separatrix E scp — 2a — 10) with the doublet n = 6. 

From a general point of view, the evolution of the en- 
ergy levels of a spectrum as a function of a given param- 
eter (here e) can be compared to the time evolution of a 
gas of particles obeying dynamical laws of the Calogero- 
Moser system type (pl|-[]64[|). This system is Hamilto- 
nian and the interaction between particles (eigenvalues) 
is strongly repulsive at short distances. This gives rise to 
avoided crossings. Nevertheless, in the presence of sym- 
metries, the different parity sectors decouple from each 
other and thus, do not interact |35|. As a result, the four 
sectors {s,s,a,a} evolve individually according to the 
Calogero-Moser dynamics but do not interact with each 
other. Thus they may cross. Such a crossing is observed 
in Fig. U in the vicinity of the point (e = 5, E = 44). 
Fig. |§| also shows the pairing of eigenvalues whose mo- 
mentum related parities are complementary. For instance 
eigenvalues of the s sector cluster with eigenvalues of the 
a or a sectors. They are indeed quasi-indistinguishable. 
The reason is that the influence of the onsite potential be- 
comes weak and thus the corresponding splittings small 
as the energy becomes sufficiently high . 

By following the quadruplet which starts at E ~ 41.43, 
we observe that it first survives an avoided crossing with 
one (|s), | a)) pair {e ~ 0.7). Then it survives again an- 
other one with two (|s), \ a)) states and starts to split into 
two pairs (|s), \a)) and (|s), \a)) after a collision with a 
quasi-quadruplet which originates from the octuplet at 
E ~ 42.3. Finally, in the vicinity of e = 7.5, the (|s), |o)) 
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FIG. 9: Evolution of |$) (|<& m ,n| 2 ) as the coupling strength 
e increases. Each picture represents the 2D Fourier space, 
pi(m) along the x axis and P2(n) along the y axis. The value 
of each component |$ m ,n| 2 is represented by its gray level 
from white (0) to black (~ 0.3) . 



pair seems to join a new (|s), |a}) pair but remains clearly 
separated from it. The fact that pairs of permutation- 
ally related eigenvalues may survive avoided crossings is 
already known and may for instance be observed in the 
case of the trimer problem ( []55| , ) . 

In order to see the progressive evolution of our initial 
state |$) as the coupling increases, we have plotted snap- 
shots of its evolution for different values of e (Fig. g). 
Each picture represents the 2D discrete Fourier space. 
The a;-axis represents the coordinate m € Z related to 
the momentum of the first pendulum pi and the y-axis, 
n € Z related to pi- The |$) state has been computed on 
a square lattice of 41 x 41 Fourier components. For clarity 
only a 21 x 21 lattice has been displayed. We observe that 
the state |$) remains almost unchanged for < e < 2. 
The excitation is well localized in the 2D Fourier space. 
The average momentum (Pi) = J2 m m\$ m . n \ 2 , which 
defines the projection of the excitation "center" on the 
x-axis, is large and makes the first pendulum "rotating" 
whereas (P2) = J2 m n n \^m,n\ 2 is small and makes the 
second pendulum " oscillating" . The quotes indicate that 
this correspondence refers to an interpretation in terms 
of the classical system. As e becomes larger than 4, the 
state starts to spread on the lattice and does not corre- 
spond to a "coherent" excitation anymore. A plot of the 
Poincare section of the classical system at e = 6 with 
an energy corresponding to the one of |<I>) has shown the 
allowed region of the Fourier space to be chaotic (except 
in the vicinity of the in-phase and out-of-phase motion 
and in a tiny region where the classical rotobreather still 
exists). This means that any small perturbation of the 
initial conditions of the classical rotobreather leads to a 
chaotic trajectory. We verify in this particular example 



that the strong overlap of the chaotic sea with the quan- 
tum state leads to large splittings of tunneling pairs (here 
(| S ),|s»and(|a>,|a»). 



5. Evolution of the splittings with e 

To provide quantitative results concerning the behav- 
ior of the different splittings involved in the expressions 
( p7| ) as e increases, we have computed them over the same 
range of values. The result is presented in Fig. [h]. It 
always concerns the state |$) obtained from the quadru- 
plet of preceeding section. Notice that for the sake of 
clarity, only three of the six possible splittings have been 
displayed on the figure. The omitted splittings behave 
very similarly according to the symmetry properties of 
the involved states. Each of the displayed splittings cor- 
responds to a given tunneling process as specified in the 
figure. After a crossing value e ~ 0.15, the graph basi- 
cally consists of two curves. One is related to the momen- 
tum or energy transfer between the pendula and another 
one is related to the reversal of the total momentum of 
the system. This means that the upper curve is related 
to the transfer of the excitation from a general point of 
view (energy or momentum). The lower curve represents 
a process (the global momentum reversal) which is not 
associated to any kind of excitation transfer between the 
pendula but only to a global modification of the system. 
It has no relation with the tunneling of the quantum ro- 
tobreather from site to site. 

Different splittings 
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Total momentum, 
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FIG. 10: Dependence of different splittings of the quadruplet 
whose energy at e — is E ~ 41.43 on e. Only three of 
them have been displayed, each being associated with a given 
tunneling process. The splittings are given in terms of the 
rescaled energy AE = 2AE. 



Note that at the anti-continuous limit, the energies of 
the states \a) (\s)) and \a) (|s)) are identical. The split- 
ting which corresponds to the energy transfer (\E a — E a \ 
in the figure) is thus zero. This explains the behavior 
of the corresponding (thin solid) curve, in logarithmic 
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scale, which decreases regularly with s. Moreover, due to 
this degeneracy, all the remaining splittings are equal and 
given by @. For the |$) state we obtain S 7 ~ 3.15 10~ 10 
which corresponds to the saturation value obtained nu- 
merically for \E S ~ E a \ and \E S — E a \ (thick solid and 
dashed lines). 

The splittings depend smoothly on e except for val- 
ues of the coupling where crossings and avoided crossings 
take place (see Fig. ||). The first crossing involves the 
quadruplet states {|s), |a), |s), |a)} and a (|s'), |a')) pair 
at e ~ 0.73. Avoided crossings occur between (|s), \s')) 
and (|a), |a')) pairs and are responsible for the peaks ap- 
pearing in the three splittings. The case of the second 
crossing (e ~ 3.30) is slightly different. This avoided 
crossing concerns the (|s), \ s')) and (\a),\a'}) pairs. The 
splitting \E S — E a \ behaves smoothly whereas the two oth- 
ers exhibit spikes due to the collision. This follows from 
the fact that \E S — E a \ does not contain any contribution 
of states interacting with the colliding pair. Note that 
the spikes observed in such crossings may have different 
forms according to the exact realization of the collision 
and the specific splitting under consideration . 

Finally, when the coupling parameter becomes larger 
than 5, the splittings corresponding to the transfer of the 
excitation between the pendula start to be of the same 
order as the mean spacing in the spectrum (1/ir here as 
shown in (j34|)). This is the signature of the disappearance 
of the quantum rotobreather since no quantum state is 
able to keep the energy on a given pendulum during a 
sufficiently long time at the considered energy level. 
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FIG. 11: Dependence of the splitting governing the energy 
tunneling on e for different quadruplet energies. The three 
curves correspond to quadruplets starting from a ground 
state and a doublet at e = 0. The respective labels are 
m = 7, 11 and 17. The corresponding energies at e — 
are E = 41.43, 113.3 and 281.2. 

We conclude this section by discussing the dependence 
of the splitting governing the energy transfer on the 
quadruplet energy (Fig. |ll|). As expected the splittings 



decrease with increasing quadruplet energy. This can be 
qualitatively understood by referring to the analytical 
expression of the splitting at a = 0. Indeed, taking into 
account the slow change of the quadruplet energy with e, 
E m ^ ~ m 2 + e 2 /(m 2 — 1), we may directly replace m by 

\J E m fi in the expression of the splitting. We thus obtain 

A(E, a = 0) ~ 4(e/2)^/r(v^E) 2 which gives a rapidly 
decreasing function of the quadruplet energy for values 
of e < E. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

In this paper, after a brief review of the essential re- 
sults of the single pendulum problem, we have used a 
two coupled pendula model to show the possibility of 
constructing a quantum state |$) whose properties are 
similar to those of a classical rotobreather. This state 
has been built in a very similar manner to those used 
in the corresponding classical system, namely, by start- 
ing from a state which mimics the behavior of a rotat- 
ing pendulum decoupled from another one at rest at the 
anti-continuous limit. Four states, each belonging to a 
different symmetry sector, are shown to be necessary to 
form it. By switching on the coupling between the pen- 
dula, we follow the resulting evolution of the quadruplet 
and monitor the different splittings. Each of them can 
be shown to be associated to a given tunneling process 
according to the involved states. 

In general, we have seen that two processes have to be 
distinguished. The first concerns the tunneling (or trans- 
fer) of the initial excitation (energy) between the pendula 
which corresponds to a tunneling of the quantum roto- 
breather |$) between neighboring sites. This tunneling 
effect is associated to the permutation symmetry present 
in the system. The second tunneling effect relies on an- 
other symmetry due to the invariance of the system to a 
reversal of its total momentum. This symmetry already 
appears in the single pendulum problem where it is re- 
sponsible for the occurrence of doublets of nearly degen- 
erated states at sufficiently high energy (above the sepa- 
ratrix). This second tunneling effect takes a time which 
is orders of magnitude larger than the time of the exci- 
tation transfer between pendula (except for very small 
coupling values). 

By progressively increasing the pendula coupling, we 
have shown that the quadruplet under investigation may 
survive crossings (or avoided crossings) with other states. 
Nevertheless, for large enough coupling, the nearly degen- 
erated states of different permutation symmetry parities 
separate sufficiently from each other, leading to splittings 
of order of the mean level spacing in the spectrum. This 
situation corresponds classically to a large chaotization 
of the phase space at the considered energy level and 
thus to a strong overlap of the quantum state with the 
chaotic layer surrounding the quasi-regular island where 
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the classical rotobreather is located. The disappearance 
of the quantum rotobreather can be interepreted by the 
fact that the corresponding state is no more able to keep 
the excitation in a given pendulum for a long time as 
compared to the typical oscillation time of the system. 



To conclude this paper, we would like to comment 
on the relation between the classical and the quantum 
(roto)breathers. An essential ingredient for the exis- 
tence of a quantum breather is the appearance of bands 
of nearly degenerated eigenstates whose bandwith re- 
main very small as the coupling increases. Nevertheless, 
the general method employed to construct the quantum 
breather by starting from the anti-continuous limit may 
be used as well for a chain of N purely harmonic os- 
cillators (with a harmonic coupling). As this chain can- 
not possess any kind of classical breathers, what prevents 
quantum breathers from existing on such a chain? To an- 
swer this question, we compute the bandwith correspond- 
ing to a local excitation of n bosons close to the anti- 
continuous limit. We obtain AE n = n{\J iS 1 + 4e — uj), lo 
being the oscillator frequency and s the coupling. This 
shows that the bandwith becomes large as the mean en- 
ergy of the band E n ~ nuj does. Moreover, the density 
of states scales like g(E) ~ E N ~ X (see the introduction) 
and the product g(E n )AE n — > oo as E — + oo . This shows 
that like the classical harmonic chain, the quantum sys- 
tem cannot support quantum breathers. 



A second comment concerns the " classical-like" behav- 
ior of the quantum breather. By defining the latter as 
a superposition of the eigenstates making up a single 
band, its time evolution is de facto restricted to these 
states. The fine structure of the band thus provides the 
only available frequencies of this evolution. As the cor- 
responding splittings are very small, none of them is re- 
lated to the classical breather frequency. These frequen- 
cies only concern the tunneling of the excitation from 
site to site. On the other hand, the individual energies of 
particles making up a classical breather, for example, os- 
cillate around their mean value at the breather frequency. 
Where does such a frequency appear in the quantum sys- 
tem? It turns out that this frequency is naturally recov- 
ered by considering quantum states which display an ex- 
citation similar to the quantum breather but which are 
not restricted to a single band. We have verified in our 
system that a coherent state parametrized by the average 
phases and momenta of the quantum rotobreather |$) ex- 
cites mainly the quadruplets separated by an energy dif- 
ference which corresponds approximately to the classical 
rotobreather frequency. The interaction between quadru- 
plets is thus responsible for the " classical- like" time be- 
havior of averaged observables (like individual energies). 
This should generalize to nonlinear systems with more 
than two degrees of freedom where we expect interactions 
between bands to play a similar role as the interaction 
between quadruplets in our system. 



APPENDIX A: COUPLED ROTORS 

In this appendix, we solve the problem of two coupled 
rotors, i.e. the limit of the two coupled pendula when 
the onsite parameter a is equal to zero. In this case, 
the system defined by the Hamiltonian (|2^) is integrable. 
Indeed, 

H CI = i (pj + pi) + e(l - cos(a;i - x 2 )) (Al) 

Let s — [x\ + x 2 )/2 and d — {x\ — X2)/2. Then the 
classical equations of motion reads 



s = 

d + £sin(2d) = 



(A2) 
(A3) 



The system decouples and consists of a free rotor for the 
center of mass motion and a pendulum for the relative 
coordinate 2d (with an "onsite" parameter 2e). 
Let us turn t o th e associated quantum problem. The 
Hamiltonian ( |Al|) gives rise to the time independent 
Schrodinger equation 



1 



-~ (C + O + e(l ~ cosOn - x 2 )) - E 



ip(xi,x 2 ) = 
(A4) 

which yields (after the change of variables defined above) 
1 



(<9 2 + dj) + e(l - cos(2d)) - E 



ip{s,d) = (A5) 



where ip(s, d) — ip(x\, x 2 ). Because of the 27r-periodicity 
of ip(xi, X2) in each of its variables we may expand it as 
a double Fourier series: 

ip(xx,x 2 ) = — ^ ' ( / , m : nexp(i(rna;i + nx 2 ))(A6) 



m,n£2 2 



s exp(z(cxs + 3d)) 



<r,8e( 



E 

er,<5e(2Z+i; 

i>(s, d) 



i>al exp(i(o-s 



3d))} (A7) 



where we have defined a = m + n and 6 = m — n and 
where the notations a, 5 G (2Z) 2 and a, S € (2Z+1) 2 indi- 
cate respectively that it, 5 have to be both even numbers 
or both odd numbers. Moreover, ^\ = V ; (o-+5)/2,(ct-5)/2 
for a, S both even (reason of the superscript ( e )) and sim- 
ilarly = ip( a +5)/2,(a-8)/2 I0r Cj 5 both odd. It follows 
that the first term in (A7) is a 7r-periodic function in 



each of the variables s or d whereas the second one is 
27r-periodic. 
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Notice further that because the potential only depends 
on d, (|A5| ) decouples to give 



d 2 - 2E S 



Ms) = o 



d\ + 2e(l - cos(2d)) - 2E d 



(A8) 

Md) = (A9) 



where ijj(s,d) = ili s (s)ipd(d) and i? = E s + Ed- 

The last equation of (A8) is a Mathieu equation with 
a characteristic value a and a parameter q defined by 



4(Ed - e) and q = —2e 



(A10) 



Taking into account the preceeding remarks concerning 
the periodicity as well as the two possible parities (odd 
or even) of the solutions, we finally get 



s 4 



and E 5 , = - 



(All) 



where a, 5 £ Z 2 have to have the same parity (even or 
odd) and where we have introduced a spinor notation 



if 
if 



A« > 
n < 



(A12) 



For the usual notation of the Mathieu equation see |J1YJ 
By defining the rescaled energy by E = 2(E—e) and usinj 
the quantum numbers m, n we finally get the followinj 
expression for the spectrum of the coupled rotors 



E m , n = - I (m + n) 2 + 




(m, n) G Z 2 
(A13) 



It follows E„ 



E. 



Thus all the levels are twofold 



degenerated except when m = —n, which corresponds 
to a total momentum of the system equal to zero. No- 
tice that, when s — > 0, the characteristic values converge 
to (m — n) 2 and we recover the free rotors spectrum 



differences 



m 



Moreover, at sufficiently high momenta 
n\ > y/e, the energies -EV^n and E n ^ rn 



become nearly degenerated, 
by means of (glj) gives 



Their splitting, computed 



AE\ 



\m— n\ 



-AE\ 



| m — n | 



z j2\\m-n\ 



(|m-n|-l)! 2 



(A14) 



The eigenfunctions corresponding to the spectrum (A13) 
may be expressed according to their symmetry sectors by 



(A15) 



where the normalization factor N a is 1 /tt for a ^ and 
1/V27T if a = and where a, 5, s and d have already 
been defined as functions of m, n, x\ and x%. C a (s) — 
cos(crs) and S a (s) = sin(crs). With this notation, the 
correspondence between the symmetry sectors and the 
values of a and S are 



\s) => (a > ; S > 0) 
|s) (a < ; S < 0) 



\a) 
IS) 



(d < ; <y > 0) 
(a > ; <5 < 0) 
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For r = 1, the correction to the leading order has to be 
computed separately. It reads —q 3 /32. This expression 
is in exact agreement with the small q expansion of the 
characteristic values. 
[69] It is interesting to compare this result with the one we 
may obtain by evaluating the number of states below the 
separatrix N s ~ 8-ya/V, a 3> 1 (see ( p7j ) in the limit 
k — > 1). The number of (nondegenerated) pairs of sym- 
metric and antisymmetric eigenstates below the separa- 
trix is thus n' c ~ TV.,/2 = 4-y/a/V which is very similar to 
n c as 8e/n ~ 1.067. As En s — E s — 2a it proves that 
expression (M) starts to be valid immediately as E > E s . 



